Event-by-Event Fluctuations in Ultrarelativistic Heavy-Ion Collisions 
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Motivated by forthcoming experiments at RHIC and LHC, we study event-by-event fluctuations in 
ultrarelativistic heavy-ion collisions in participant nucleon as well as thermal models. The calculated 
physical observables, including multiplicity, kaon to pion ratios, and transverse momenta agree well 
with recent NA49 data at the SPS, and indicate that such studies do not yet reveal the presence of 
new physics. Finally, we present a simple model of how a first order phase transition can be signaled 
by very large fluctuations. 
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Central ultrarelativistic collisions at RHIC and LHC 
are expected to produce at least ~ 10 4 particles, and 
thus present one with the remarkable opportunity to an- 
alyze, on an event-by-event basis, fluctuations in physical 
observables such particle multiplicities, transverse mo- 
menta, correlations and ratios. Analysis of single events 
with large statistics can reveal very different physics 
than studying averages over a large statistical sample of 
events. The use of Hanbury Brown-Twiss correlations to 
extract the system geometry is a familiar application of 
event- by-event fluctuations in nuclear collisions [[l], and 
elsewhere, e.g, in sonoluminesence Q. The power of 
this tool has been strikingly illustrated in study of in- 
terference between Bose-Einstein condensates in trapped 
atomic systems ||. Recently NA49 has presented a pro- 
totypical event-by-event analysis of fluctuations in cen- 
tral Pb+Pb collisions at 158 GeV per nucleon at the 
SPS, which produce more than a thousand particles per 
event ||. 

Studying event-by-event fluctuations in ultrarelativis- 
tic heavy ion collisions to extract new physics was pro- 
posed in a series of papers by the authors and co- 
workers [5-7], in which the analysis of transverse energy 
fluctuations in central collisions |§ was used to extract 
evidence within the binary collision picture for color, or 
cross-section, fluctuations. More recent theoretical pa- 
pers have focussed on different aspects of these fluctu- 
ations, such as searching for evidence for thcrmaliza- 
tion and critical fluctuations at the QCD phase 

transition |T|Q. 

In order to be able to extract new physics associated 
with fluctuations, it is necessary to understand the role 
of expected statistical fluctuations. Our aim here is to 
study the sources of these fluctuations in collisions. As 
we shall see, the current NA49 data can be essentially 
understood on the basis of straightforward statistical ar- 
guments. Expected sources of fluctuations include im- 
pact parameter fluctuations, fluctuations in the number 
of primary collisions, and the results of such collisions, 
nuclear deformations S , effects of rescattering of secon- 



daries, and qcd color fluctuations. Since fluctuations in 
collisions are sensitive to the amount of rescattering of 
secondaries taking place, we discuss in detail two limit- 
ing cases, the participant or "wounded nucleon model" 
(WNM), in which one assumes that particle production 
occurs in the individual participant nuclcons and rescat- 
tering of secondaries is ignored, and the thermal limit 
in which scatterings bring the system into local ther- 
mal equilibrium. Whether rescatterings increase relative 
fluctuations through greater production of multiplicity, 
transverse momenta, etc., or decrease fluctuations by in- 
volving a greater number of degrees of freedom, is not 
immediately obvious. Indeed VENUS simulations |l3| ] 
showed that rescattering had negligible effects on trans- 
verse energy fluctuations. As we shall see, both models 
give similar results for multiplicity fluctuations. In the 
wounded nucleon model fluctuations arise mainly from 
multiplicity fluctuations for each participant and from 
impact parameter fluctuations. Limited acceptance also 
influences the observed fluctuations. We calculate in de- 
tail statistical fluctuations in multiplicity, K/7T ratios, 
and transverse momentum. Finally, we show in a sim- 
ple model how first order phase transitions are capable 
of producing very significant fluctuations. 

Multiplicity fluctuations: Let us first calculate fluc- 
tuations in the participant model, which appears to de- 
scribe well physics at SPS energies B. In this picture 
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where N p is the number of participants and rii is the 
number of particles produced in the acceptance by par- 
ticipant i. In the absence of correlations between N p and 
n, the average multiplicity is (N) — (N p )(n). For exam- 
ple, NA49 measures charged particles in the rapidity re- 
gion 4 < y < 5.5 and finds (N) ~ 270 for central Pb+Pb 
collisions. Finite impact parameters (6^S3.5 fm) as well 
as surface diffuseness reduce the number of participants 
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from the total number of nucleons 2A to (N p ) ~ 350 es- 
timated from Glauber theory; thus (n) ~ 0.77. Squaring 
Eq. (Q) assuming (niUj) = (n,i)(nj) for i ^ j, we find the 
multiplicity fluctuations 



u) N = u n + (n)ui Np , 
where in general we write 

p{y) = (y 2 ) - (y) 2 = (y)v 



(2) 



(3) 



for any stochastic variable y. 

A major source of multiplicity fluctuations per partic- 
ipant, w„, is the limited acceptance. While each partic- 
ipant produces v charged particles, only a smaller frac- 
tion / = (n) j (v) are accepted. Without carrying out a 
detailed analysis of the acceptance, one can make a sim- 
ple statistical estimate assuming that the particles are 
accepted randomly, in which case n is binomially dis- 
tributed with a(n) = vf(l — /) for fixed v. Including 
fluctuations in v we obtain, similarly to Eq. (J2j), 



u„ = 1 - / + fuJ v 



(4) 



In NN collisions at SPS energies, the charged particle 
multiplicity is ~ 7.3 and oj v ~ 1.9 113; thus {v) ~ 3.7 
and / ~ 0.21 for the NA49 acceptance. Consequently, 
we find from Eq. (^) that u n ~ 1.2. 

As a consequence of nuclear correlations, which 
strongly reduce density fluctuations in the colliding nu- 
clei, the fluctuations w/v p (6) in N p are very small for fixed 
impact parameter b 0j. Almost all nucleons in the nu- 
clear overlap volume collide and participate. [By con- 
trast, the fluctuations in the number of binary collisions 
is non-negligible.] Cross section fluctuations play a small 
role in the WNM |Q . Fluctuations in the number of par- 
ticipants can arise when the target nucleus is deformed, 
since the orientations of the deformation axes vary from 
event to event jl5| |. The fluctuations, ujn p : m the num- 
ber of participants are dominated by the varying impact 
parameters selected by the experiment. In the NA49 ex- 
periment, for example, the zero degree calorimeter selects 
the 5% most central collisions, corresponding to impact 
parameters smaller than a centrality cut on impact pa- 
rameter, b c ~ 3.5 fm. We have 



"n p (N p ) = \ (° d 2 bN p (b) 2 - (N p ) 

7r0 c JO 



(5) 



where (N p ) = (l/nb 2 c ) d 2 bN p (b). The number of 
participants for a given centrality, calculated in [pj| , 
can be approximated by N p (b) ~ N p (0)(l — b/2R) for 
< fe<3.5 fm; thus 



N p {0) ( b, 
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For NA49 Pb+Pb collisions with N p {0) ~ 400 and 
(b c /2R) 2 ~ 5% we find lon p — 1-1. Impact parameter 



fluctuations are thus important even for the centrality 
trigger of NA49. Varying the centrality cut or b c to con- 
trol such impact parameter fluctuations (JsJ) should enable 
one to extract better any more interesting intrinsic fluctu- 
ations. Recent WA98 analyses confirm that fluctuations 
in photons and pions grow approximately linearly with 
the centrality cut The Gaussian multiplicity distri- 
bution found in central collisions changes for minimum 
bias to a plateau-like distribution Q . 

Calculating u>m for the NA49 parameters, we find from 
Eq. (§), cu N ~ 1.2 + (0.77)(1.1) = 2.0, in good agreement 
with experiment, which measures a multiplicity distri 
bution oc cxp[-(7V - {N)) 2 /2{N)oj 
order 2.01 @. 

Let us now consider, in the opposite limit of consid- 
erable rescattering, fluctuations in thermal models. In 
a gas in equilibrium, the mean number of particles per 
bosonic mode n a is given by 



c.l'/n 
J N ] 



where is of 



(exp(£ Q /T)-l) 



(7) 



with fluctuations 



uj na = 1 + (n > . (8) 
The total fluctuation in the multiplicity, N ~ ^2 a n a , is 



, , be 

N 



If the modes are taken to be momentum states, the result- 
ing fluctuations are lu^ e = £(2)/£(3) = 1.37 for massless 
particles, while for pions at temperature T = 150 MeV 



LU% E = 1.11 |18| 



Resonances add to fluctuations in the thermal limit 
whereas they are implicitly included in the WNM fluc- 
tuations. In high energy nuclear collisions, resonance de- 
cays such as p — > 2n, uj — > 37r, etc., lead to half or more of 
the pion multiplicity. Only a small fraction r ~ 10% pro- 
duce two charged particles in a thermal hadron gas |l9| l 
or in RQMD (^0|. Including such resonance fluctuations 
in the BE fluctuations gives, similarly to Eq. (||), 



1 - r 



+ (1 + r)u>* E 



(10) 



With r ~ 0.1 we obtain lu b ! e+r ~ 1.3. If not all of the 
decay particles fall into the NA49 acceptance the fluctu- 
ations from resonances will be reduced. In Jl^] the esti- 
mated effect of resonances is about twice ours: lom ~ 1.5, 
not including impact parameter fluctuations. 

Fluctuations in the effective collision volume add a fur- 
ther term (N)a(V)/(V) 2 to to% E+R . Assuming that the 
volume scales with the number of participants, uoy / {V) ~ 
w Np /(N p ), we find from Eq. (§) that u N = u EE+R + 
{n)u>M p — 2.1, again consistent with the NA49 data. Be- 
cause of the similarity between the magnitudes of the 
thermal and WNM multiplicity fluctuations, the present 
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measurements cannot distinguish between these two lim- 
iting pictures. 

Kaon/pion ratio: To second order in the fluctuations 
of the numbers of K and tt, we have 

The corresponding fluctuations in (K/tt) are given by 

(Kn)-(K)(ir) 



D 2 = ^K/tt 



(K/tt) (K) (tt) 



- 2- 



(12) 



The fluctuations in the kaon to pion ratio is dominated 
by the fluctuations in the number of kaons alone. The 
third term in Eq. ( |l2|) includes correlations between the 
number of pions and kaons. It contains a negative part 
from volume fluctuations, which removes the volume fluc- 
tuations in lux and u T since such fluctuations cancel 
in any ratio. In the NA49 data Q the average ra- 
tio of charged kaons to charged pions is (K/tt) = 0.18 
and (tt) ~ 220. Excluding volume fluctuations, we take 
lok — lo-k — 1.2 — 1.3 as discussed above. The first two 
terms in Eq. ( |l2| ) then yield D ~ 0.20 in good agreement 
with preliminary measurements D = 0.23 [Q. Thus at 
this stage the data gives no evidence for correlated pro- 
duction of K and 7r, as described by the final term in 
Eq. (12), besides volume fluctuations. The similar fluc- 
tuations in mixed event analyses P> m ixed — 0.208 B con- 
firm this conclusion. 

Transverse momentum fluctuations: The total 
transverse momentum per event 



N 



Ft = ^Pt,i , 



(13) 



is very similar to the transverse energy, for which fluc- 
tuations have been studied extensively [|]||. The mean 
transverse momentum and inverse slopes of distributions 
generally increase with centrality or multiplicity. Assum- 
ing that a = dlog((p t ) N)/d log N is small, as is the case 
for pions [pl| , the average transverse momentum per par- 
ticle for given multiplicity N is to leading order 



(pt)N = (pt)(l + a(N-(N))/(N)). 



(14) 



where (pt) is the average over all events of the single par- 
ticle transverse momentum. With this parametrization, 
the average total transverse momentum per particle in 
an event obeys (P t /N) = (p t ). When the transverse 
momentum is approximately exponentially distributed 
with inverse slope T in a given event, {pt,i) — 2T, and 
(j{pt,i) = 2T 2 = (p t ) 2 /2. This latter fluctuation is in 
principle dependent on multiplicity, but as a higher or- 
der effect, we ignore it in the following. 

The total transverse momentum per particle in an 
event has fluctuations 



(N)*(P t /N) = a(p t ) + a 2 ( Pt ) 2 Lu N 

^E(^m-^> 2 )>- ( 15 ) 



The three terms on the right are respectively: 

i) The individual fluctuations <j{pt,i) — (Pti) ~ (Pt) 2 , 
the main term. In the NA49 data, (p t ) = 377 MeV 
and (N) = 270. From Eq. @ we thus obtain 

{aiPt/N) 1 ' 2 /^) ~ 1/V2p) = 4 - 3 %: which accounts 
for most of the experimentally measured fluctuation 
4.65% The data contains no indication of intrinsic 
temperature fluctuations in the collisions. 

ii) Effects of correlations between p t and N, which are 
suppressed with respect to the first term by a factor ~ a 2 . 
In NA49 the multiplicity of charged particles is mainly 
that of pions for which T ~ (pt)/2 increases little com- 
pared with pp collisions, and a ~ 0.05 — 0.1. Thus, these 
correlations are small for the NA49 data. However, for 
kaons and protons, a can be an order of magnitude larger 
as their distributions are strongly affected by the flow ob- 
served in central collisions pl| . 

iii) Correlations between transverse momenta of differ- 
ent particles in the same event. In the WNM the mo- 
menta of particles originating from the same participant 
are correlated. In Lund string fragmentation, for exam- 
ple, a quark-antiquark pair is produced with the same 
Pt but in opposite direction. The average number of 
pairs of hadrons from the same participant is (n(n — 1)) 
and therefore the latter term in Eq. ( |l5| ) becomes 
({n(n - l))(n))((p t:i p t j^i) - (p t ) 2 )- To a good approx- 
imation, n is Poisson distributed, i.e., {n(n — l))/{n) 

[11 . equal to 0.77 for the NA49 acceptance, so that this 
latter term becomes ~ ((pt,iPt,j=a) — (pt) 2 )- The mo- 
mentum correlation between two particles from the same 
participant is expected to be a small fraction of aOpt^). 

To quantify the effect of rescatterings, Ref. M sug- 
gested studying the differences in (N)a(P t /N) and cr{pt) 
via the quantity 



$(p t ) ~ y/(N)<T(P t /N) - yJ<T(pt,i) ■ 



(16) 



As we see from Eq. (|15|), in the applicable limit that the 
second and third terms are small, 



*(pt) 



1 



(a 2 (p t ) 2 uj N + {{p t ,iPt,jjti) ~ (Pt) 2 )) ■ 

(17) 



In the Fritiof model, based on the WNM with no rescat- 
terings between secondaries, one finds $(pt) — 4.5 MeV. 
In the thermal limit the correlations in Eq. ( fHf ) should 
vanish for classical particles but the interference of iden- 
tical particles (HBT correlations) contribute to these cor- 
relations by ~ 6.5 MeV Jl(J and slightly reduced by res- 
onances. The NA49 experimental value, Q(pt) = 5 MeV 
(corrected for two-track resolution) seems to favor the 
thermal limit Q. Note however that with a ~ 0.05 — 0.1, 
the second term on the right side of Eq. (|7h alone leads 
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to $ ~ 1 — 4 MeV, i.e., the same order of magnitude. 
If ((pt,iPt,jjti) — (Pt) 2 ) is not positive, then one cannot a 
priori rule out that the smallness of $(pt) does not arise 
from a cancellation of this term with a 2 (p t ) 2 LUN, rather 
than from thermalization. 

First order phase transitions can lead to rather 
large fluctuations in physical quantities. Thus, detec- 
tion of enhanced fluctuations, beyond the elementary sta- 
tistical ones considered to this point, could signal the 
presence of such a transition. For example, matter un- 
dergoing a transition from chirally symmetric to broken 
chiral symmetry could, when expanding, supercool and 
form droplets, resulting in large multiplicity versus ra- 
pidity fluctuations Q . Let us imagine that Np droplets 
fall into the acceptance, each producing n particles, i.e., 
(N) = (Nu)(n). The corresponding multiplicity fluctua- 
tion can be computed analogously to Eq. (0) 



uj n = uj n + {n)uj ND . 



(18) 



As in Eq. (0), we expect u> n ~ 1. However, unlike the 
case of participant fluctuations, the second term in JlS| ) 
can lead to huge multiplicity fluctuations when only a 
few droplets fall into the acceptance; in such a case, 
(n) is large and lon d of order unity. The fluctuations 
from droplets depends on the total number of droplets, 
the spread in rapidity of particles from a droplet, Sy ~ 
■^/T/mt, as well as the experimental acceptance in rapid- 
ity, Ay. When Sy Ay and the droplets are binomially 
distributed in rapidity, uin d ~ 1 — Ay/y tot , which can be 
a significant fraction of unity. 

In the extreme case where none or only one droplet 
falls into the acceptance with equal probability, we have 
luno = 1/2 and (n) — 2(A). The resulting fluctuation 
is lon ~ (N), which is more than two orders of magni- 
tude larger than the expected value of order unity as cur- 
rently measured in NA49. This simple example clearly 
demonstrates the importance of event-by-event fluctu- 
ations accompanying phase transitions, and illustrates 
how monitoring such fluctuations versus centrality be- 
comes a promising signal, in the upcoming RHIC exper- 
iments, for the onset of a transition. 
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